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THE MIRRORS MODEL : MACROSCOPIC DIFFUSION WITHOUT NOISE OR CHAOS 


YANN CHIFFAUDEL AND RAPHAEL LEFEVERE 


Abstract. Before stating our main result, we first clarify through classical examples the 
status of the laws of macroscopic physics as laws of large numbers. We next consider 
the mirrors model in a finite d-dimensional domain and connected to particles reservoirs 
at fixed chemical potentials. The dynamics is purely deterministic and non-ergodic but 
takes place in a random environment. We study the macroscopic current of particles in 
the stationary regime. We show first that when the size of the system goes to infinity, 
the behaviour of the stationary current of particles is governed by the proportion of orbits 
crossing the system. This allows to formulate a necessary and sufficient condition on the 
distribution of the set of orbits that ensures the validity of Tick’s law. Using this approach, 
we show that Tick’s law relating the stationary macroscopic current of particles to the 
concentration difference holds in three dimensions and above. The negative correlations 
between crossing orbits play a key role in the argument. 


1. Macroscopic laws as laws of large numbers 

Take a macroscopic box A = [0, L]'^ that contain N freely moving distinguishable par¬ 
ticles and fixed obstacles of arbitrary shapes. N should he thought to be of the order of 
magnitude of the Avogadro number : 6 x 10^^. 

A first experiment is performed on this box. The N particles are initially located in a cube 
A' c A of side length L' < L, see Figure The evolution of the density of the cloud 



Figure 1. In the first experiment, the cloud of particles is initially concen¬ 
trated in a volume A'. 

of particles is monitored through a beam of light that crosses the system. We call this 
density p : A x [0, (X)[— )• p(x, t) G M+. The initial state is described by p(x, 0) = 1a'(x)/|A'|. 
The empirical fact that is observed at the macroscopic level is that the density evolves 
according to the laws of diffusion: 

( dtp{x,t) = KAp{x,t) 

\ • Vp(x, t) = 0, X G 9A (1.1) 

i P(x,0) = po(x) := 1a'(x)/|A'|, 

where rix is the vector normal to the boundary of the box dA at x and k is a strictly positive 
constant. How can we explain this phenomenon from the motion of the individual atoms ? 
For each macroscopic coordinate x = (xi,..., Xd) G A, we define a microscopic coordinate 
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q = r/eAT where en = The motion of the particles is entirely determined by the 

law of Newtonian mechanics. When a particle makes a collision with one of the fixed 
obstacles or the boundaries of the boxes, its velocity is modified according to the laws of 
specular reflection. We assume that each particle starts with a speed equal to 1. Since this 
property is preserved by the dynamics, the microscopic motion of a given particle (with 
label i G {!,... ,N}) is described by a map t —)> (qj(e^^f), pi(e^^t)) G [0,L/€n]^ x 
where is the unit sphere in d dimensions. The coordinates (q*, pj) are the microscopic 
positions and velocities of the i-th particle. The microscopic time-scale is The scaling 
of the time variable is a priori arbitrary but is fixed here by the fact that the solution of the 
diffusion equation p{x, t) is invariant under the transformation (x, t) —)> (Ax, X^t), A > 0. 

In the absence of any other information, we assume that the initial positions and veloc¬ 
ities of the particles are independent and identically uniformly distributed, with density 


/(qi,Pi) 


pi 

|^/||grf-l| lA'(e7Vq^)lg>^-i(p») 


( 1 . 2 ) 


for every i = 1,... ,N. If the position of each particle is chosen independently of the 
others with that density, the number of particles in a microscopic volume of size of order 
1 follows a Poisson distribution with finite mean as W —)> oo. We denote by P the law of 
probability of the initial positions and velocities of particles. No information about the 
spatial location or the shape of the obstacles in A is known either. We denote by Q the 
probability distribution on those degrees of freedom. It is chosen such that in the limit 
N ^ oo, the number of obstacles in a microscopic volume of size 1 follows a distribution 
with a finite mean and such that, almost surely, the dynamics of moving particles is well- 
defined at all time. The dynamical system defined in this way is an instance of the random 
Lorentz gas. 

We define the empirical density of particles : 


1 ^ 

PN{x,t) = 


(1.3) 


j=i 


The density p^ contain all possible information about the density. Indeed, it is easy to see 
that 

PNiV,t) := / dx ly(x)pAr(x,t) 

dRd 

gives the proportion of particles that belong to any 1/ c A at time t. It is straightforward 
to see that the following statement holds : if {(qj(0), pi(0)) : z = 1,..., is a collection 
of i.i.d variables with marginals given by (|1.2p then for any bounded function h and any 
<5 > 0 : 

lim P[| {pn{ 0), h) - (po, /i) | > <5] = 0, (1.4) 

N^oo 

dx /i(x)p(x) and we use the notation pNit) ■= 


where po is given by (1.1), {h,g) = 

Pn{-, t)- The goal is to show the 

Conjecture 1.1. There exists a natural^distribution Q such that for any t > 0, any bounded 
function h and any 6 > 0 : 


lim ] 

N^OO 


X Q[| {PN{t),h) - {p{t), /i) I > d] = 0. 


where p{t) := p{-,t) is the solution of ([L J| for some k > 0. 


(1.5) 


1 


By natural we mean as uniform as possible over the locations and shapes of obstacles. 
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The law of ordinary diffusion is therefore understood as a law of large numbers: as N 
becomes very large, the probability that the empirical density pN{t) differs significantly 
from the solution of the diffusion equation goes to zero. 

It is natural to consider first a simpler version of the problem in which the randomness 
of the obstacles is removed, i.e. Q is taken to be a Dirac distribution 6c on a special 
configuration of obstacles giving rise to a chaotic dynamics. This is exactly the result of 
Bunimovich and Sinai ffSl. They consider the 2D case in which obstacles are disks located 
at the vertices of a regular lattice such that the induced billiard dynamics has a finite 
horizon]^ Their result implies that for any t > 0, any bounded function h and any <5 > 0 : 

lim P X (Jed {pNit), h) - {p(t), /i) I > (5] = 0. (1.6) 

N^oo 

Let us sketch how this statement is obtained. Let /i : —)• M a bounded function. First, 

one computes : 


1 iV 

{pN{t),h) = ^ dx V 5 (eArqi(e^^f) - x)/i(x) 

1 ^ 


(1.7) 

( 1 . 8 ) 


i=i 

Thus, because the initial positions of the particles are identically distributed : 

E[{pN{t),h)] = E[h{eNCii{effH))]. 

Next, the theorem 2 of [l3]| implies]^ that 

lim E[/i(eArqi(eT^t))] = / p{x,t)h{x,t) dx 

N^OO ./nd 


(1.9) 


where p(x, t) is the solution of (1.1). To derive (1.9), one has to rely on the strong chaotic 


properties of the billard system under study. 

Next, since {qj(f) : 1 < j < are independent, the variance of (/Oiv(t), h) is 


Yav[{pN{t),h)] = ^Var[/i(eArqi(e 


^*))l = 0(1) 


since h is bounded. Thus, Chebychev inequality allows us to conclude the proof of Con¬ 
jecture |1.1| in the case where Q = 6c- One should note that the proof is made of two 


steps of very different levels of complexity. The first step is basically given by ( |1.9 ) and 
this is where the whole difficulty is located. The second step is a concentration result of 
the random variable (/>Ar(t), h) around the expected value E[{pN{t),h)]. This part is trivial 
because when Q is replaced by 6c, the positions and velocities of the particles remain 
independent for all time t. The fact that it is so trivial is probably the reason why it is 
hard to find a reference where this step is mentioned or even alluded to. It is however 
essential and when Q / 6c, the statistical independence of the motions of particles is 
lost. Controlling the correlations between them to ensure the concentration of {pN{t), h) 
around its mean does require some work. We will see below that this issue arises in the 
mirrors model and how it can be dealt with. 

A second experiment may be performed on the box. At the two sides of the cube A 
perpendicular to ei = (1,0...,0), particles reservoirs maintain constant values of the 
local densities of particles pc and p^, respectively on the left and right side, see Figure 2. 


Q 

For instance, the center of each (sufficiently large) disk is located at a vertex of a triangular lattice. 

To be more precise Bunimovich and Sinai consider the case L = oo but their method should apply directly 
to the finite L case 
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Figure 2. At the two sides of the cube A, particles reservoirs maintain 
constant values of the local densities of particles pi and p^. 

A device records the net flux of mass crossing a section of A perpendicular to ei per unit 
time. This quantity is denoted by j{x, t) when the section contains the point (x, 0,..., 0) 
for X G [0, L]. After some transient time proportional to it is observed that the instan¬ 
taneous current of particles takes the stationary value : 

3s{x) = ‘^{PL- Pr). (1.10) 

With respect to the first experiment, the coupling to external reservoirs introduce an 
additional probabilistic element. The law of the reservoirs and the initial conditions of 
the particles inside the system is denoted by P. One can introduce an empirical current of 
particles per unit time J 7 v(t)|^ Again, the goal is to show the following conjecture : 

Conjecture 1.2. There exists a natural distribution Q such that for any bounded continuous 
function h and any J > 0 : 

lim lim P x Q[| (jN{t),h\ — {js,h) | > <5] = 0. (1-11) 

N^oot^oc \ / 


While this has never been done explicitly, one should expect that with the choice Q = (5c 
the result follows from the methods of flS]] . In Dll , the authors show that in a low density 
regime limit, the expectation of the stationary current (with respect to P x Q ) converges 
to js- The statement corresponding to Conjecture |1.2| together with the exponential con¬ 
vergence to the stationary current has been obtained in the case of a discrete space-time 
dynamics in DSl. We now outline how the problem may be tackled in the context of the 
mirrors model. 


2. The mirrors model 

The mirrors model was introduced by Ruijgrok and Cohen |0 as a lattice version of 
the random Lorentz gas or the Ehrenfest wind-tree model. The latter encompasses a large 
class of models in which obstacles do not induce a chaotic behaviour of the trajectories 
of the particles. A fundamental question which remains open regarding those models 
is whether a non-chaotic deterministic dynamics may give rise to a macroscopic diffusive 
behaviour. In the mirrors model, particles travel on the edges of the cubic lattice generated 
by “Mirrors” are located at the vertices of the lattice and deflect the motion of an 

^This quantity will be our main object of study in the next section and will be given a precise definition 
there. 
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incoming particle in a new direction, see Figure for an illustration in the 2D version of 
the model. The precise general definition is given below. 



Figure 3. In the mirrors model on I?, green particles travel on the edge 
of the lattice and are deflected by mirrors located at the vertices. The green 
particles do not interact with each other. 


A quick look at the structure of the orbits reveals its total lack of ergodicity. Indeed, in 
Figure]^ a sample of orbits in a finite box with periodic boundary conditions in the vertical 
direction and reflecting boundary conditions in the horizontal direction is pictured with 
different colors. For almost any configuration of the mirrors, no orbit is able to visit the 
entire phase space. 

A perhaps even more striking fact is that, in any dimension, the motion of a particle 
in an environment of randomly orientated mirrors is not a gaussian diffusion []4l] . More 
precisely this means that ( 1.9| ) (where the expectation is taken with respect to P and Q) 
does not hold. 

Our goal is to show that in spite of these unpromising properties, the mirrors model 
does exhibit normal macroscopic conductive properties when d > 3 in the sense that the 
analogue of ( 1.11 ) holds. It turns out that quite weak conditions on the statistics of orbits 
are sufficient to ensure the validity of Fick’s law at the macroscopic level. It is therefore 
not necessary that orbits behave as a Gaussian diffusion to ensure the validity of Fick’s law. 
Thus, the normal macroscopic laws of diffusion apply to a much wider class of dynamical 
systems than generally expected. 

The dynamics of the mirrors model is reversible in the usual sense of the word in the 
context of Hamiltonian dynamics. Namely, under the reversal of the velocities of all parti¬ 
cles at a given time t > 0, the dynamics brings the system of particles to its initial condition 
at time 0 (with reversed velocities), see ( |2.2 ). This reversibility property of the dynamics 
will allow us to show that when the system is large, the number of orbits travelling from 
one side of the system to the other one basically determines the value of the current in the 
stationary state. This will allow to formulate a condition on the distribution of orbits that 
is both sufficient and necessary for the validity of Fick’s law. 


^ In spite of this, it has been observed numerically @ that in two dimensions, the mean-square displace¬ 
ment of a given particle is linear in time, allowing the definition of a microscopic diffusion coefficient. This is 
of course a much weaker property than the property ra- In particular we will see that one can not define 
a macroscopic diffusion coefficient in 2D. 
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We recall now briefly the set-up of the original mirrors model. Particles travel on the 
edges of I? with unit speed. Mirrors are located at some vertices of the lattice and take 
two possible angular orientations : {f, When a particle hits a mirror, it gets deflected 
according to the laws of specular reflection, see Figure for sample trajectories of parti¬ 
cles. It is convenient to think that every particle starts at time zero with a given velocity 
at a vertex of the lattice Q that is obtained by taking the middle point of every edge of 
I?. As all particles move with unit velocity, one can simply observe the evolution of the 
system at discrete times t G N. At those times, the particles will be always located at one 
of the vertices of the new lattice Q with a well-defined velocity. In general, the orientation 
of the mirrors is picked randomly. It is obvious that the motion of a single particle can not 
be described as a Markov process. When a particle hits a mirror for the second time, no 
matter how far back in the past the first visit occurred, its reflection is strongly affected 
by the way its was reflected at the first visit. For instance in Figure the two orientations 
of the mirrors are picked at random, and in that case, at the second visit the reflection is 
always deterministic. 

We come now to a more general definition of the dynamics in d dimensions. We denote 
by z = (zi,..., Zrf) a generic element of Z'^. As for I?, we consider the set of midpoints 
of edges of an hypercuhe of Z^ of side N and with periodic conditions in all but the 
first direction. We call this set Q. It may he described as follows : Q = U^=i where 
Lj = {z + : 0 < 2:1 < A^ — 1, ( 22 , ■ ■ ■,Zd) ^ (Z/A^Z)'^”^}. Let (ei ,... ,6^) the canoni¬ 

cal basis of M'^, the space of possible velocities isV = {±^,...,±^} and the phase space 
of the dynamics is 

= {(q> p) : q G Q, p G "P S. t. if q G Lj then p = iy}- 

We denote a generic point of Af by (q, p). The set of points in M whose spatial coordinate 
belongs to the boundaries of the system is B = U P+, with 

B- = {x = e M : ci= {qi,... ,qd) e Li,qi = ^} 

B+ = {x = (q,p) G At : q = (gi,... G Li,gi = A - 

See Figure]^ 



Figure 4. The spatial component of the phase space M and of the sets P_ 
and P+ in 2D. 
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For each z G we define the action of a “mirror” on the velocity of an incoming 
particle by 7r(z; •) which is a bijection of V into itself. It satisfies the following conditions : 

( 7r(z; —7r(z; p)) = —p, Vz G Z^,Vp G P 

\ vr(0,Z2,...,Zrf;-f) = f (2.1) 

[ 7r(iV, = (Z2, ...,Zd)e (ZINZ)<^-^ 

The dynamics is defined on M in the following way. For any (q, p) G : 


^(q> p) = (q + P + 7r(q + p; p), 7r(q + p; p)). 


It is easy to check that the map F is a bijection on J\4. The two last conditions in ( |2.ip are 
just saying that when particles hit the boundaries they are reflected backwards. We define 
an operator R ■. M ^ M. which reverses the velocities by i?(q, p) = (q, —p) and we see 
that the first of the conditions (|2.1|) ensures that the map F is reversible, 




= RFR. 


( 2 . 2 ) 


We define the orbit of a point x G M., Ox = {y G M : 3t > 0, F*(x) = y} and its 
period, T{x) = inf{t > 0 : F*(x) = x}. From the fact that F is bijective, one infers that for 
every x G M., Ox is a loop : T{x) < \M\ and that orbits are non-intersecting •. \iy ^ Ox, 
then Oxi^Oy = it). A given orbit is also non-self-intersecting : if ?/ G Ox and y ^ x then 

F{y) / F{x). 

As we are interested in the transport of particles, we define occupation variables cr(q, p; t) 
{0,1} that record the absence or presence of a particle at position q with velocity p at time 
t G N. When connecting the system to external particles reservoirs, we obtain the follow¬ 
ing evolution rule : given — 1), we define cr(-; t) for all t G N* by 

( a{F~^{x);t — 1) if x ^ U 
(T(x;t) = < (T“(f —1) if X G 

1 1) if xe B+. 

The families of random variables : x G i?_, t G N} and {cr+(t) : x G F+, t G N} 

consist of independent Bernoulli variables with respective parameters and /7+. If one 
chooses {cr(x;0) : x G M} to be a collection of independent random variables, then it 
is easy to see by induction that at any t > 0, {a{x-,t) : x G M} is a collection of i.i.d 
Bernoulli random variables. To simplify a bit the discussion, we choose an homogeneous 
initial distribution, i.e. all Bernoulli random variables have a common parameter pj. The 
distribution of the collection {a{x;t) : x G M} becomes stationary after a finite time. 
More precisely, for any t > \M.\, we have the following equality in law : 

( ai if Oxr\B = % 

a{x,t) = < CT_ if F“**(x) G 

[ a+ if F-**(x)gF+ 

where t* = inf{t : F“*(x) G B} and a± and aj are Bernoulli random variables of parame¬ 
ter p± and pj. 

Proceeding as in [|8]], it is possible to show that when the size of the system goes to 
infinity, the stationary current converges in probability to the proportion of crossing orbits 
times the chemical potentials difference. We define the average current of particles that 
crosses the hyperplane Q* = {q G Q : = / + |}, Z G {1,..., — 2} during a diffusive 

time interval ; 

.j 

= Y. ^(a;;s)A(x,Z) 

s=i+l x&M 


(2.3) 
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where A{x,l) = 2(p-ei)lqgg!,with x = (q, p). Thus A{x,l) takes the value +1 (resp. —1) 
if X crosses the slice from left to right (resp. from right to left). We denote by Af± the 
numbers of crossings from B± to i?=p induced by F, i.e. Af± = I'S'il where S± is given by 

S± = {xe B±:3s> 0,V0 < j < s, F^{x) ^ B±,F\x) G B^}. 


5 _ 

4 _ 

3 _ 

2 _ 

ft 

V 

'' < 

. ' 

t y 

f y 

^ V > 

* y 

' < 
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’ 4 

f y 

V t 

‘V 



Figure 5. N = 6. Crossing orbits are coloured in purple, internal loops in 
black and non-crossing orbits are coloured in green. The travel direction 
given by the arrows is arbitrary. Each edge of the crossing orbits will be 
used twice in a given orbit; once in each direction. For this configuration 
of mirrors M = 2. 


One notes that A/+ = A/l. Indeed, since every orbit is closed, it must contain as many 
left-to-right than right-to-left crossings. Thus, we set M = A/+ = M-. Proceeding as in 
Il8l , we get that for any t > 0 Moreover, for every <5 > 0, 

any t > \M\ and / G {1,..., — 2}, 


M 

j^{p-- P+) 


> 5 


< 2exp(-<52w 


d+l \ 


(2.4) 


We take now random configurations of reflectors {7r(z; •) : z G lA}. The law of the 
reflectors is denoted by Q. The map F becomes now a random map. 

The model satisfies Pick’s law if and only if there exists some k > 0 (the conductivity) 
such that V(5 > 0, 


lim lim P X Q[| A^J((, t) — k(p_ 
N^oo t—>-oo 


P+)\ > <5] = 0. 


(2.5) 


As in DBl, it is easy to infer from (2.4) that the following theorem holds. 


Theorem 2.1. Sufficient and Necessary Condition for Pick’s law .• (2.51 holds if and 
only if there exists k > 0 such that for any <5 > 0, 


lim Q 
N^oo 


]yd-2 


> 6 


= 0 . 


( 2 . 6 ) 


^This relation implies that the average current flows in the “right” direction and that when p_ / p+, the 
average current in the stationary state is different from 0 if and only if Af / 0. 
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We see that the central object to study is the distribution of the number of crossing 
orbits M. The expectation of this quantity is related to the probability that one orbit 
crosses the system, while the variance is given in terms of the joint probability that orbits 
with two different starting points cross the system. Indeed by periodicity, we have, using 
the notations O = (( 5 ,0,..., 0), and S = S- : 


E 


■ Af ■ 

J\fd-2 


^ E[ 1 , 65 ] = Nq[0 G 5] 

xeB- 


(2.7) 


and 

( 2 . 8 ) 

with 

6{x,y) = Q[x G G S'] - Q[x G S]Q[y G 5]. (2.9) 

Thus if the two following 


Var 


Nd-2 


Crossing Conditions, are satisfied : 

(1) There exist k > 0 such that the RHS of ( [2.7J ) converges to k as N 

(2) The RHS of ( [2.8| goes to zero as N —)■ 00 . 


■ CX). 


then Pick’s law ( |2.5p holds in the stationary state. We note first that when d = 2, (2.6) 
can not hold, whatever the distribution Q is. To see this, we adapt an argument found in 
[71 . Indeed, the spatial part of each crossing orbit crosses any “vertical” section an odd 
number of times. On the other hand, the spatial part of any non-crossing orbit must cross 
any vertical section an even number of times, see Figure Thus, N and Af must have 
the same parity. This implies that there can not exist k > 0 such that (2.6) holds when 
d = 2. The origin of this issue lies in the strong correlations between crossing orbits that 
are present in two dimensions. 

We turn now to the higher dimensional case d > 3 equipped with some natural and 
spatially homogeneous distribution Q. Now observe that ifQ[7r(z; |) = -^] > 0 for some 
i = 1,... ,d then an orbit starting from O will encounter this type of reflecting mirror after 
an exponential number of steps and therefore Q[0 G 5] < for some c > 0. This, in 


turn, implies that limAr_).oo A^Q[0 G 5] = 0 and that Pick’s law can not hold. Thus from 
now on, we consider maps such that 7r(z;|) f 


zi < 


N and such that the 

conditions (|2.1|) are satisfied. We call the set of such maps 11. We take Q such that the 


if 0 < 


collection of maps 

{ 7 r(z;.) : 0 < < iV, {z 2 , ■ • •, G {ZlNZf-^} 

is independent and that each map is uniformly distributed over IT. We note first that if 
the law of an orbit with respect to Q was similar to the law of a simple random walk, 
then there would be a k > 0 such that lim 7 v-s.oo iVQ[0 G S] = k, this follows from the 
gambler’s ruin argument. Similarly, if the orbits were independent objects, then the RHS 


of (2.8) would go to zero because the only non-zero term would be the one with x = O 
and Q[0 G 5] ~ k/N. We also note that the average stationary current is identified as 
the difference between chemical potentials times the probability that a particle crosses 
the system, an idea that was put forward in [S], in the context of chaotic systems. The 
law Q of the mirrors induces a law on the set of orbits which is a priori very far from the 
distribution of independent simple random walks. The set of orbits is a very interesting 
lattice object in itself which features some (self-)avoiding properties as we mentioned 
above. 
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Fortunately, what is needed to ensure the validity of ( |2.6P is much less than the full joint 
distribution of the orbits. Thanks to (2^1 and ( |2.8P , one only has to analyze the marginal 
of a path starting on the boundary and also the joint probability of two such paths. The 
distribution of a path starting at O (i.e. on the boundary) is similar to the one of a “true” 
self-avoiding random walk DJI but defined on Q rather than on and with further con¬ 
straints. The diffusive behaviour of those walks for d > 3 has been conjectured in [l2Tl . 
see also the rigorous results of IfTOll . It can be expected that as the dimensionality of the 
system increases, the effect of the revisits of an orbit to the same mirror decreases. In a 
process where the mirrors are flipped randomly after being used (i.e memory effects are 
killed), we computed that in d = 3 the crossing probability is ~ 3/2N. Numerical simula¬ 
tions in d = 3 show that this number is indeed a good approximation. The log log plot of 
the crossing probability Q[0 G S] is given in Figure |^for N up to 400. The corresponding 
conductivity is k = 1.535 ± 0.005. As the conductivity measured in simulations is slightly 
higher than 3/2, it indicates that recollisions tend to push forward the orbit. 



Figure 6 . Q(0 g S) for N from 5 to 420. The 95% confidence interval is 
about half the size of a dot. 

We must show now that YIxgb^ 6(0, x) —> 0 as > oo. We know that this sum is 
positive because it is a variance, and thus it is enough to get an upper bound on the 
sum. We will use a numerical analysis to show that 6(0, x) < 0 for all but a finite (i.e. 
independent of N] number x G i?_. 

But before doing that and to get a better picture of the origin of the correlations d(0, x), 
we split them in two parts. Given an orbit O, we denote by 7 ( 0 ) the set of edges of Z'^ 
used by O. For each z G Zjy := {z G Z*^ ; 1 < zi < — 1, (z 2 ,... ,Zd) G (Z/A^Z)*^}, let also 

6 z(C>) be the half of the number of times that the orbit O visits the vertex z. Two crossing 
orbits O and O' are incompatible if 7 ( 0 ) n 'j(O') / 0 and compatible otherwise. The law 
of a given crossing orbit is : 

bJO) 

WO) = n n niTjfT-c 

zeZjv i=i ^ 

If Q((9, O') is the joint probability of two orbits O and O' then Q(C>, O') = 0 when O and 
O' are incompatible. If they are compatible, the joint law of two crossing orbits is given 
by 

bz{0)-F6z(0') 

q(o,o')= n n 

zeZiv i=i 


2(d-i) + l- 


( 2 . 11 ) 
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In particular, if they do not share any mirrors, then Q(C>, O') 
properties, starting from (2.9), we obtain that for x G i?_. 


Q(0)Q(0'). From those 


J(0,x) = (Q(Oo,Ox)-Q(Oo)Q(0.)) 

00 , 0 :. 

- ^ Q(Oo)Q(Ox). (2.12) 

0o,0. 


Both sums run over orbits that cross the box Q. The first sum runs over compatible orbits 
such that j(Oo) and j(Ox) share a vertex of Z'^. The second (prime) sum runs over 
incompatible orbits Oq, Ox- 

Thus, from ( 2.12^ , we see that correlations 6{0,x) are created from two opposite ori¬ 
gins, corresponding to each of the two sums in ( 2.12 ). If two orbits O and O' share some 
mirrors, then it is easy to see from (2.11) that Q(C>, O') > Q(C>)Q(C>'). This in turn im¬ 
plies that the first sum is strictly positive. This is a “cooperative” effect, the orbits help 
each other crossing the system. The second sum corresponds to a jamming effect : an 
orbit starting from O and crossing the system occupies a certain number of horizontal 
edges. Because distinct orbits can not share the same edges, the occupied edges are no 
more available for an orbit starting from x G B-, this creates negative correlations. 

Numerical simulations in d = 3 show that the latter effect dominates. For all but a few 
points, the correlations <5(0, x) for x ^ O are not only small but negative within confidence 
intervals, see Figure [7 The only exceptions are points ((1/2,1,0), ^), ((1/2,0,!),^), 
((1/2, N — 1, 0), aM ((1/2,0, iV — 1), which give clearly positive correlations. How¬ 
ever, we checked that for N = 70, ((1/2, y, 0), ^)) = —1.360 x 10“°^ ± 1.47 x 


10“°^, i.e. it is negative with a margin of more than Ou. ((1/2,0, z), ^)) must 

be equal by symmetry. Increasing values of N do not modify this behaviour. In partic¬ 
ular, the number of points with positive correlations do not increase. Since we know 
already that Q[0 G S'] ~ k/N, as N —^ cx), we conclude with the same margin that 

YIx&B- <5(0, x) < k/N — )• 0, as iV —)• cx). 



Figure 7. 5(0, x) for x = ((1/2, y, 0), ^). N = 70. We draw the 95% confi¬ 
dence interval. 


We expect the same behaviour in d > 3. A rigorous proof that the crossing conditions 
introduced above are satisfied seems to be within reach in the present model. Moreover, 
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it is possible to draw a general conclusion from the above discussion. If one is really inter¬ 
ested in deriving macroscopic laws from microscopic dynamics, many detailed properties 
of the latter are irrelevant. Only weaker properties than chaoticity ergodicity or Gaussian 
behaviour of the orbits are required. In the present context, the minimal properties neces¬ 
sary to obtain Pick’s law are encapsulated in the crossing conditions. It is of course natural 
to seek similar weak conditions in different contexts as for instance in the problem of the 
derivation of Fourier’s law. 
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